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' Abstract. We consider an 50(4) Euler rigid body with two 'inertia momenta' coinciding. 

, We study it from the point of view of bihamiltonian geometry. We show how to algebraically 

pL^ ' integrate it by means of the method of separation of variables. 
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> ' 

in 

O ■ 1 Introduction 

\^ ' It is fair to say that the problem of the free rigid body in (i-e., the SO (A) Euler or Euler- 

■ Manakov rigid body or top0) can be still considered as an interesting problem in the theory of 

i Separation of Variables (SoV) for Hamilton-Jacobi equations. 

Dh' In genious methods were devised to solve the very classical problem of the Euler-Manakov 

^ ■ top, (see, e.g., [H [H [3l O [191 ED ESI [27j, but this hst - as well as the bibhography of the present 

"j^ , paper - is by far incomplete), in particular to characterise the solutions by means of suitable 

techniques coming from algebraic geometry. 

In this paper we will consider, from the point of view of bihamiltonian geometry, a degenerate 
case of Euler top that can be called rotationally symmetric. Namely, we consider the case in 
which two of the four independent parameters that enter the construction (that, after Manakov's 
' construction, can be termed generalised principal inertia moments) coincide. This is a particular 

case of the general 50(4) system, but it possesses some interesting and peculiar features that 
deserve, in the author's opinion, to be spelled out. 

It should be noticed that in [26] the SoV problem was implicitly considered in the study of SoV 
for elliptic Gaudin models, and basically solved, with a considerable amount of hard computa- 
tions and ingenuity. Furthermore, in [20j . interesting classes of solutions to the Hamilton-Jacobi 
equations associated with the SO {4) Euler-Manakov systems (as well as other alike systems) 
were recently obtained. In the present paper we will follow a more direct and geometrical way to 
approach it, with the aim of giving a simple and explicit setting for the Separation of Variables 
problem of this symmetric SO (4) Euler top. 

The framework we will use to study this system is the so-called bihamiltonian setting for 
Separation of Variables for Gel'fand-Zakharevich |12j systems that has been introduced a few 
years ago (see, e.g., [HO [23], and references quoted therein) and formalised in [10]; we will 
sketch the content of this method in Section [2j 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 



The full collection is available at http://www.emis.de/journals/SIGMA/kuznetsovlSmll 
^Actually, this model was studied, in the XIX century, by Schottky and Kotter. 
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After that, in Section [3] we will briefly resume those aspects of the SO{4:) Euler-Manakov 
system that are relevant in our analysis, and in the core of the paper (Section UD we will consider 
the rotationally symmetric case and solve, applying the recipes described in Section [2l the SoV 
problem for this Hamiltonian system. 

Separation of Variables, both in classical and quantum systems, was a research arena in which 
Vadim Kuznetsov obtained notable results. In particular, his papers [131 [TSl [HI [16] were deeply 
influential on the study of the SoV problem (and more general integrability) from the point of 
view of bihamiltonian geometry, carried out by the author of the present paper in collaboration 
with F. Magri, M. Pedroni, and G. Tondo. We had the opportunity to discuss with him these 
and related subjects many times, and appreciate his scientific, as well as human talents. With 
deep sorrow I dedicate this work to his memory. 

2 A bihamiltonian set-up for SoV 

In this Section we will briefly review a scheme for solving the SoV problem in the Hamilton 
Jacobi equations, based on properties of bihamiltonian manifolds. We will discuss those features 
that are relevant for the case at hand, referring to [10] and to [9] for a comprehensive theoretical 
presentation, as well as for a wider list of references. 

Let {M, Pi, P2) be a bihamiltonian manifold, that is, a manifold endowed with a pair of 
compatible Poisson brackets {■,-}pi, i = 1,2, or, equivalently, with two compatible Poisson 
bivectors (or operators) Pi, P2, related to the brackets by the well-known formulas 

{f,g}p, = {df,P,dg) yf,geC^{M), i = l,2. 

We consider a Gel'fand-Zakharevich bihamiltonian system with one bihamiltonian chain. That 
is, we consider the datum, on the bihamiltonian manifold (M, Pi, P2), of an 'anchored' Lenard- 
Magri chain of length n > 0, 

PidHo = 0, PidHi = P2dHi_i, i = l,...,n, P2dHn = 0, (2.1) 

that is, a family of bihamiltonian vector fields originating from a Casimir function of one Poisson 
operator and ending in a Casimir function of the other Poisson operator. We may suppose that 
p additional Casimir functions common to the two structures Ci, . . . ,Cp be also presenlH, and 
we require that the system be complete, i.e., that n = dimM — 1—p independent vector fields 
fill in the chain (f2T]) . 

Such a system provides (families of) Liouville integrable systems as follows. One considers a 
(generic) symplectic leaf S of one of the two Poisson operators, say. Pi. 

5 is a submanifold of M defined by fixing the values of all Casimir functions Hq, Ci, . . . ,Cp 
of Pi. Any of the vector fields of the chain (j2.ip . say Xi = PidHi{= P2dHQ)), restricts 
to S, is still Hamiltonian with respect to (the restriction of) the chosen Poisson operator Pi 
(which becomes an ordinary symplectic operator on S), and, thanks to the basic property of 
Lenard Magri chains, comes equipped with the right number of involutive integrals, namely 
the (restriction to S of) the other Hamiltonian functions of the vector fields Xi, i = 2, . . . ,N 
of CT . 

What is lost, in general, in this procedure, is the bihamiltonian structure of the equations: 
indeed, the second Poisson operator does not restrict to S, since it does not leave the function Hq, 
which is a Casimir of Pq invariant, i.e., Hamitonian vector fields generated by P2 do not leave 
the sub manifolds S invariant. 

However, as it is shown in [10], the manifold M may be provided with another suitable 
bihamiltonian structure. That is, along with Pi, a "new" second Poisson operator Q can be 



^In the 5*0(4) Euler top, indeed, we will have one of such Casimir functions. 
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defined deforming P2; Q can be restricted to the symplectic leaves S and has those properties 
needed for integrating the vector fields of the chain (|2.1|) . 

For this new structure to exist, conditions are to be fulfilled. Namely, one has to find a vector 
field Z, defined on M, such that: 

i) Z is a symmetry of Pi, transversal to the submanifolds = const, and leaving the 
common Casimir functions Ci invariant: 

Liez(Pi) = 0, UeziHo) = l, Liez(C,) = 0, i = l,...,p. 

ii) It deforms of the second Poisson tensor Pi as follows: 

Liez(P2) =Y A Z, for some vector field Y. 

Indeed, under these conditions it holds that the bivector Q := P2 — Xi A Z (where Xi is the 
first vector field of the Lenard-Magri chain (12. ip ) satisfies: 

1. Q is a Poisson structure on M, compatible with Pi which shares with Pi all the Casimirs, 
and hence, together with Pi induces a bihamiltonian structure on the symplectic leaves S. 

2. The Hamiltonians Hi do not form anymore Lenard-Magri sequences w.r.t. this new Poisson 
pair {Pi,Q), but still are in involution also w.r.t. the deformed (or new) structure Q = P2 — 
XiAZ. 

These two properties are very important for our purposes; indeed, from the first one it follows 
[m [I7j that the pair {Pi,Q) defines, on each symplectic leaf S, a special set of co-ordinates, 
called Darboux-Nijenhuis co-ordinates, associated with the eigenvalues of a torsionless 'recursion' 
operator. 

From the second one, it follows [10] that the Hamiltonians Hi are separable in these DN co- 
ordinates, that is, that the Hamilton- J acobi equation associated with any of the Hamiltonians Hi 
is separable in these co-ordinates. Hence, the first and basic step of the bihamiltonian recipe 
for SoV of Gel'fand-Zakharevich systems of type (12. Ih essentially boils down to find/guess this 
vector field Z, which will be referred to as the transversal vector field. 

For the reader's convenience, as well as to provide the necessary background to the calcu- 
lations presented in Section HI we discuss more in details the construction of DN co-ordinates. 
A preliminary remark is in order. As the two structures Pi and Q share the same symplectic 
leaves, we will generically use the same letters Pi, Q also for their natural restrictions to the 
symplectic leaves to avoid cumbersome notations. 

As already noticed, on any (generic) symplectic leaf of Pi this operator is symplectic and 
thus invertible; hence, the compositions 

N := QP~^ : TS TS, and N* := P^^Q : T*S T*S 

are well defined. In the literature, A'^ is called a Nijenhuis or recursion or hereditary operator; 
in our setting, its adjoint operator, A^* will play a more visible role. 

Being the ratio of two antisymmetric operators, N* has up to m = (1/2) dim 5 distinct 
eigenvalues Ai, . . . , A^- Under the assumption that the number of these distinct eigenvalues 
be exactly m, it follows, basically as a consequence of the compatibility condition between Pi 
and Q, that, for each Aj there is a pair of canonical (w.r.t. Pi) co-ordinates fi, gi that generates 
the eigenspaces of A*, that is, it holds 

N*dfi = Xidfi; N*dg„ {fi, gj}p, = Sij. (2.2) 

These co-ordinates are called Darboux-Nijenhuis co-ordinates associated with the pair Pi, Q. 
Notice that the 'eigenvalue' relations written above imply the Poisson bracket relations {fi,gj}Q 
= 6ijXi w.r.t. the second Poisson structure. Fortunately enough, in generic cases, there is no 
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need to integrate (all the) two-dimensional distributions Ker(A^* — Aj) to actually find DN 
co-ordinates, thanks to the following results. 

It should be noticed that, since A'^* depends on the point of S, its eigenvalues are functions 
defined on S. In particular, non-constant eigenvalues of N* do provide DN co-ordinates. This 
means that, if dXi ^ 0, then one can choose the function Aj as the co-ordinate function /j 

of a . 

To find the missing canonical co-ordinate associated with the eigenvalue Aj one can try to 
use a recipe, discussed in [10] that might be called method of deformation of the Hamiltonians, 
and goes as follows: 

Consider the sum pi of all the eigenvalues of N* , and the Hamiltonian vector field Y = 
—Pidpi- Then collect the Hamiltonians filling the Lenard-Magri recursion relations (|2.ip in the 
Gel'fand-Zakharevich polynomial 

H{X) := X^Ho + A"-^i7i + • • • + XHn^i + (2.3) 

and deform it repeatedly along the vector field Y, that is, consider the polynomials 

H'{X) = LieyF(A), H"{X) = Liey//'(A), .... 

If, for some n > it happens that the polynomial H^'^~^'^\X) identically vanishes in A, (while 
H^^~^^\X) is not identically vanishing), then evaluating the rational function 

^ ^ at A = A,- 



(A) 

provides us with a DN co-ordinate gi conjugate with fi = Aj. 



Remarks: 

1) When the eigenvalues Aj are no more functionally independent, no general algorithms/recipes 
are known (at least to the author of the present paper) to find Darboux-Nijenhuis co-ordinates. 
However, when the Nijenhuis operator has a constant eigenvalue, say, A„, a case that happens 
for the symmetric /SO (4) top herewith considered, one can proceed as follows. One can at first 
apply the method briefly recalled above to construct a set of (n — 1) pairs of Darboux-Nijenhuis 
co-ordinates, associated with the n — 1 co-ordinates {Aa}JJzj defined by the Nijenhuis tensor. 
Then one has to consider the distribution associated with the constant eigenvalue A„, and try 
to integrate it (that is, to find the functions /„ and gn of equation (j2.2p ) by ad hoc methods. 

2) It should be stressed that, with respect to the pair Pi,Q, the Hamiltonians do not fill in 
a standard Lenard-Magri chain, but rather a generalised one (see [HJ [lOl [18] for further details) 
of the form 

QdHi = P^dHij^i + PidHi, i = 1, . . . , n, 

where pi are (up to signs) the elementary symmetric polynomials associated with the eigenval- 
ues Aj. However, this invariance relation is sufficient to ensure separability of the Hamilton- 
Jacobi equations. 

3) Contrary to other methods for integrating of Hamilton equations, and notably the method 
of Lax pairs, the bihamiltonian setting herewith briefly sketched provides, on general grounds, 
somewhat poor information on the Jacobi separation relations ~ that is, the relations tying pairs 
of separation co-ordinates with the Hamiltonians Hq , Hi , H2 , ■ ■ ■ , Hn and the common Casimir 
functions C = (Ci, . . . , Cp). 

However, the functional form of the separation relations can be sometimes ascertained from 
bihamiltonian geometry. Indeed, if the second Lie derivative of the GZ polynomial (j2.3p with 
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respect to the transversal vector field Z vanishes, then these relations will be affine functions of 
the Hamiltonians and the Casimirs, that is, they will be given by expressions of the form 

Fl{Xi,Ci)Hi + --- + F^{X^,^^)Hn + G{Xi,^^;C) = 0, i = l,...,m. 

Separation relations of this kind are often referred to as generalised Stackel separation relations. 

3 The Euler— Manakov model 

In this section we will briefly review the basic features of the 5*0 (4) Euler-Manakov top. 

The phase space is the (dual of) the Lie Algebra 50 (4), identifieql with 4x4 antisymmetric 
matrices 

4 

M = ^ mij{Eij - Eji), 
i<j=i 

where Eij is the elementary matrix with 1 at the (i,j)-th place. 

This six dimensional manifold is naturally endowed with the Lie Poisson structure, that, in 
the natural variables m = { mi^2T ^1,3,1^1,4, ^2,3, ^2,ii '^3,4 } is represented by the matrix 






-"12,3 


-7712,4 


"T-1,3 


7711,4 





m-2,3 





-"13,4 


-7711,2 





7711,4 


7712,4 


"13,4 








-7711,2 


-"^1,3 


-"11,3 


"11,2 








-"^-3,4 


7712,4 


-mi,4 





"11,2 


"T-3,4 





-"12,3 





-"11,4 


"il,3 


-7772,4 


m2,3 






The Hamiltonian is the quadratic function 

i<j=i 

where the coefficients Uij can be written as 

ttij = jf + j|, with {i, j, I, k} a permutation of {1, 2, 3, 4}. 

The Hamilton equations of motion (that is, the Euler equations for the 50(4) rigid body), are 
quadratic equations in the variables rriij that depend parametrically on the coefficients jf . For 
instance, 

^77112 = Jl^(77li,377l2,3 + "ll,4"l2,4) - -^2^ ("11,37772,3 + "ll,4"l2,4) 

and so on and so forth. 

Complete Liouville integrability of the model is ensured following well known facts. 

^We are actually identifying so(4) and its dual. 
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1. The rank of the so (4) Lie Poisson structure is 4; its Casimir functions are 



^o = E 



C = mi oma 4 + jni 47722,3 - r77i 37712,4- 



i<j 



2. The second independent non-trivial constant of the motion for is provided by the 
quadratic function 



i<j = l 



■-'I "^fc) 



with {i, J, Z, A;} a permutation of {1, 2, 3, 4}. 



The Hamiltonian vector field X associated with admits a Lax representation with para- 
meter [31 [19]; indeed, if one considers the matrix J := diag(Ji, J2, J3, J4), and forms the matrix 

L(A) = AJ2 + M, (3.2) 

the Euler equations are equivalent to the Lax equations 

|l(A) = [L(A),S(A)], 

where B{\) = VL + and is the 'matrix of angular velocities', i.e., defined by M via 
M = SVL + ^]. 

As it is well known, the integrals of the motion (as well as the Casimirs of the Lie Poisson 
structure) are collected in the characteristic polynomial of L{\). In particular, if one uses the 
product p\ as 'eigenvalue parameter', one gets 



Det(L(A) - pAl) = \\Pa{p)) + y{p\H^) + pHi + H2) + C\ 



where Pa{p) = H {^1 ~ P)-, Hi = —2Hs, H2 = Kg, and Hq, C are the two Casimirs of the Lie 



Poisson structure (the second one being the Pfaffian of M. 

As it was discovered in [6l [7], and independently in [22], the Euler-Manakov equations of 
motion admit a bihamiltonian formulation that can be described as follows. The matrix 
defines a deformed commutator on the Lie algebra so(4) as: 



[Ml, Mi]j2 := [Ml J2 M2, M2 Ml] = Mi]^M2 - M2J^Mi. 



(3.3) 



The Lie Poisson structure associated with [•,-]j2 provides a second Hamiltonian structure P2 
on so (4). Compatibility with the standard one is assured by the method of augmented transla- 
tions, i.e., by the fact that [•, ■]x = [•, •]j2 — A[-, •] is a one-parameter family of commutators, that 
is, the Jacobi identity holds identically in A. The 6x6 matrix representing the second Poisson 
structure in the phase space variables m = {7i7jj}i<j=i,...,4 can be easily found to be 







Jl "l2,3 



Jl 7772,4 

-■h'^rni^s 

-J2^7ni,4 





-Jl 7772,3 


W.l,2J3^ 



-"71,4^3^ 



-Jl 7772,4 J2 mi' 



-Jl 7773,4 




7771, 2 J4^ 

•mi^zJi^ 



-777i,2J3' 






-7772,4J3^ 



J2 7771,4 


-777i,2J4^ 



7772,3-'4^ 





777i,4J3^ 
-"7l,3J4^ 

7772,4J3^ 
-7772,3 J4^ 
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In particular, the so(4) Euler system is a Hamitonian vector field also w.r.t. the Poisson ope- 
rator P2, defined by the deformed commutator (|3.3|) . with 'second' Hamiltonian the function 
— {1/2)Hq = —(1/2) ''^ij- Moreover, a direct computation ensures the following: 

i<j 

Proposition 1. The characteristic polynomial (|3.3p V{X,p) = Det(L(A) — {p\)l) of the Lax 
matrix L(A) is a Casimir of the Poisson pencil Pp = P2 — pPi, that is, 

P2{dV{\p)) = pPi{dV{\,p)) identically in p, A. 

In other words, the Pfaffian C is a common Casimir of the two structures, while the three 
functions Hq, Hi, H2 satisfy the GZ recurrence relations 

Pid{Ho) = 0, P2d{Ho) = Pid{Hi), P2d{Hi) = Pid{H2), P2d{H2) = Q. 



For the sequel of the paper, the following well known considerations will be useful. 
The Lie algebra 50 (4) is isomorphic to the direct sum 50 (3) ©so (3) ; a linear change of variables 
that explicitly realises this isomorphism is the following: 

111 

a^i = - m.3,4), yi = -j={mi;i+m2,A), zi = —(mi,4 - 771,2,3), 

3^2 = ^("7-1,2 + ^3,4), y2 = ;^('Tll,3 - 777,2,4), ^^2 = ("7,1,4 + ^,2,3) ■ 

In particular, the variables {xi,yi, Zi}i=i^2 satisfy, with respect to the standard Lie Poisson 
structure Pi the so (3) commutation relations: 

{xi,yi}pj_ = V2zi, {xi,Zi}p^ = -V2yi, {yi,Zi}p^ = V2xi, i = 1,2, 

while Poisson brackets involving co-ordinates from different so (3) subalgebras vanish, e.g. 
{xi, Z2}p^ = and so on and so forth. 

Under this co-ordinate change, the Euler Hamiltonian (j3.ip acquires the form 

Hg = 2p,iXiX2 + 2p3yiy2 + 2^2^12:2 + Pi {vi' + Vi^ + X2 + xx + zx Z2) , (3.4) 
where the new constants pi are related with the Jj's as follows: 

J\ = -pi + pi- P3- P2, J2^ = P'i- Pi + Pi+ P2, 
J3^ = pi - P3 + P4 + P2, Ji^ = -P2 + Pi + P3 + Pi- 
One can notice that the Hamiltonian (|3.4p is the sum of a multiple of the Casimir function Hq 
of the standard Lie-Poisson structure, and the classical analogue 



HxYZ = 2(^4x2x1 + P3yiy2 + P2Z1Z2) 



of the Hamiltonian of a (2-site) XY Z Heisenberg model. 

We will, in the remainder of the present paper, consider a special case, namely the one that 
goes under the name of XX Z model, that is we will study the case p^x = p^. From the point 
of view of the Euler rigid body in so (4), this is tantamount to consider a rigid body with two 
principal inertia moments (J| and J|) that are equal. 



8 



G. Falqui 



4 The symmetric (or XXZ) Euler systems 

In the case ^4 = ^3 , the non-trivial part of the Hamiltonian reads 
Hxxz = 2 f^3ix2Xi + yiy2) + 2^22^1 ^2- 

This fact suggests to choose hnear co-ordinates in g = 50 (3) © so (3) adapted to the sym- 
metries of Hxxz] the choice we will follow in the sequel will be to consider the sixtuple 
{ui,vi, zi,U2,V2, Z2} related with the standard so(3) ©so(3) co-ordinates {xk,yk, Zk}i=i,2 by 



Uk = Xk + iyk, 



Vk = Xk- lyk, 



= 1,2. 



In these co-ordinates, the characteristic polynomial of the Lax matrix L[\) associated with the 
problem (that is, the one given in (j3.2p . with fJ-i = IJ.3) has the expression 

\^{Pa{p)) + xHp^Ho + pHi + W2) + \cl 

where 

Ho = uivi + V2U2 + z\' + Z2^ , C2 = U2V2 + 2:2^ - UiVi - Zi^, 
= -2/U3('U2t;i + V2U1) - AP2Z1Z2 - 2p{Hq, 

H2 = PlHo + A1I1P2Z1Z2 + 'ip?,{pi+ lJi2){v2Ul+U2Vl) + plizi' + Z2^ - VlUi - V2U2) 

- 2filizi - Z2tf. (4.1) 
The explicit expressions of the Poisson tensors (in the new co-ordinates) are, respectively. 









-A2 



with 



2zi 
-2zi 

Ui -Vi 



-Ui 





(4.2) 



and P2 = A*i -Pl + A, with Pi stiU given by (|iT2]) . and 

-U2 



A = fl2 



2z2 


U2 



* 






Ui 





-Vl 


-V2 





-2zi 

















+ PS 







V2 2{zi - Z2) 

-m 


* 



- zi> 


-U2 





V2 


Vl 








Ul 





-Vl 








(We indicate with a * the lower diagonal part of these antisymmetric tensors.) 

The Hamiltonian vector field Xi (i.e., up to a numeric factor, the Euler-Manakov 50(4) 
field in the rotationally symmetric case), generated under Pi by the Hamiltonian TCi is explicitly 
given, in these new co-ordinates, by 



d d d 

Xl = A{p2UlZ2 - H3U2Zl)^ AilJ'2VlZ2 - P3V2Zi)t:^ 4(/i2U2Zl - P3UiZ2)^ 

OUi OVi OU2 



d 

+ ^P2V2Zl - fl3ViZ2)-^ h 2fl3{u2Vl - U1V2) 

OV2 



d d 
dzi dz2 



(4.3) 



We will study the SoV problem for this Hamiltonian vector field, within the scheme of SoV for 
GZ systems of [10], resumed in Section [2l Namely we have to: 
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1. Find a suitable transversal vector field Z; 

2. Consider the Poisson operators Pi,Q := P2 — Xi A Z, as well as their restrictions to the 
generic symplectic leaves; 

3. Find the DN co-ordinates associated with these restrictions; 

4. Find the Jacobi separation relations linking pairs of DN co-ordinates and the Hamiltonians. 
One can check that a suitable transversal vector field for the problem is given by 

Z = —— + — — 

2ui dvi 2u2 dv2 

Namely, Z is a symmetry of Pi, that is, Liez{Pi) = 0. Moreover, 

LiezCHo) = 1, Liez(C2) = 0, 
and one indeed can check that the bivector 



Q = P2-XiAZ 



turns out to be a (generically rank 4) Poisson operator compatible with Pi that admits TCi 
and C2 as Casimir functions, that is, shares with the Poisson tensor Pi, associated with the 
standard Lie algebra structure of so (4) = 50 (3) ©50 (3) the same symplectic leaves. 

Remark. For further use, we notice that a direct computation shows that the second Lie 
derivative of the characteristic polynomial Det(L(A) — {pX)l) w.r.t. the transversal vector field Z 
vanishes as well. 



Liez (Liez(Det(L(A) - (pA)l)) = 0, 



(4.4) 



that is, the condition for Stackel separability is fulfilled. 

From the theoretical framework recalled in Section [2] we know that the symplectic leaves of Pi 
are four dimensional manifolds endowed with a pair of compatible Poisson structure, i.e., the 
restrictions of Pi and of Q. 

These four dimensional symplectic leaves S are obtained fixing the values {'Ho,C2} of these 
common Casimir function. Furthermore, thanks to the explicit expressions given in the first 
line of ()4.ip . in (open sets of) these symplectic leaves one can use, as co-ordinates, the four 
parameters 

U = {ui,Zl,U2,Z2}, 

since one can express the co-ordinates vi,V2 as follows: 



Vl = - 
2 



1 no-C2-2zi^ 



V2 = — 

2 



1 Wo + C2 - 2^2 



U2 



The restrictions P and Q of the Poisson structures Pi and Q to the leaf S are represented by 
4x4 matrices that have, in these co-ordinates, the explicit expressions 







-ui 


* U2 





Q 



-{fJ^3U2 + fJ-lUl) IJ.2UI - fJ'3U2 

IJ,2U2 - fisUl 

* IJ.1U2 + fisui 
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The transpose Nijenhuis operator N* = P is given by 

/W3 hMi +^2 

Ui Ml 

/is h/n 

Ui 

Ui Ui 

IJ-2 fJ'3 h/il 

/i2 /X3 h/Ui 

^2 U2 



^^2 

/is /^2 



(4.5) 



Its eigenvalues are 

Al = /XI + /i2, 



A2 = /il - /t2 + /X3 ^ • 

U2 Ui ' 



From the general theory, we know that A2 is one of the DN co-ordinates we are looking for 
(which we need to complement with its conjugate co-ordinate ^2)) while we need to find both 
canonical co-ordinates relative to the constant eigenvalue Ai. 

The problem of finding the canonical co-ordinate conjugated to the non-constant eigenvalue 
A2 can be dealt with the idea of deforming the Hamiltonian polynomial. Indeed we consider the 
sum pi of the eigevalues of N, = 2/ii -|- fJ,3{ui/u2 + U2/U1), and the vector field Y = —Pdpi; 
it is given by 



Y = /.3G(u) 



d d 
dzi dz2 



where the function G{u) is given by 

^/ N U2 Ui 

G(u) = , 

Ul U2 

and is connected with the eigenvalue A2 by 

2 



G(u) 



4 + 



A2 - Ml + /^2 
/^3 



(4.6) 



(4.7) 



Now we iteratively apply the vector field Y to the polynomial containing the relevant Hamilto- 
nians, that is 

n{p) = p^nQ + pni + n2. 

By means of direct computation one can check that Liey('H) factors as 
LieH^) = '^^(^-^^-^^^ G(u)L(u) 



Ul U2 



with 



L(u) = /is [z2Ui^ + ZiU2^) - P.2U1U2 {Zi + Z2) . 

Since, quite obviously Liey(ni U2) = 0, and, thanks to (j4.7p . 

Liey(G(u)) ={G(u),A2}p = 0, 

we are left, for the computation of the second Lie derivative of 7i, with the computation of 
Liey(L(u)). It gives 



Liey(L(u)) = /i3G(u)(uiti2-F(u)), 
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where 



F{u) 



-2^2 



A2 - (m + ^2)- 



U1U2 

So, the third Lie derivative of TC w.r.t. Y vanishes, and so the function 



6 



LieYin) 



Liey (Liey (W)) 



1 



P=A2 



L(u) 



fxsui U2 VG'(u)F(u) 



is the DN co-ordinate conjugated to A2 we were looking for. 

The two functions F(u) and G{u) wih play a role in the last task we will deal with, that is, 
the determination of the Jacobi separation relations. 

Next we turn to consider the problem of finding DN co-ordinates associated with the constant 
eigenvalue Ai = ^1 -|- ^2 of the tensor A^* of equation (|4.5p . Being this a constant eigenvalue, 
we have to find 'by hands' the associated DN co-ordinates. Fortunately enough the expression 
of the operator N* — Ai is given by 



-/^2 + 



U2 



Ui 







/^3 







-/^2 + 



U2 



Ui 



^2^2 
Ui 



U2 



M3 



M2 











-/^2 + 



/^3^1 
U2 



fJ'3U2 
Ui 



Its kernel can be easily found to be generated by the two 1-forms 

ai = dz2 - dzi, a2 = (usUi - ^2U2)dui + {^3^2 - ^Ji2'Ui)du2 
that integrate, respectively, to the functions 

Ci 



1 



1 



Z2 - 21, 



2^3 ^il 



H2U1U2 + -/U3 U2 



A direct computation shows that {Ci^^iIp = — 2^i, and so the DN co-ordinate conjugated to Ci 
is ^1 = — (1/2) log ^1. Thus, we have found, on the generic common symplectic leaf S of Pi 
and Q, the desired set of DN co-ordinates (Ci, = —(1/2) log 9i, A2, ^2)- It can be noticed that, 
along with the two Casimirs Ci, C2, they provide a set of co-ordinates in a Zariski open subset 
of the phase space so (4) ~ so (3) © so (3). 

What we are left with the determination of the Jacobi separation relations, namely we have 
to seek for two relations of the form 



^'i(H;Ci,ei) = 0, <I>2(H;A2,6) = 0, (4.8) 

linking pairs of DN co-ordinates and the conserved quantities H = {C2, "^i, "^2} of (j4.1|) . 

Owing to the different ways the separation co-ordinates were found, and, especially, the fact 
that one of the separating momenta is an additional constant of the motion, we expect that $1 
and $2 have different functional dependence on their variables; so, instead of trying to use 
the spectral curve relations we directly seek for the Jacobi relations (j4.8p . by means of explicit 
calculations. Since the characteristic condition for Stackel separability is verified in our case 
(see equation (j4.4p ). we can look for Jacobi relations that are affine functions in the Casimirs 
and the Hamiltonians. 

Let us first consider ^i(H; Ci, Ci) = 0; i^i this respect, one can notice that (^i is an additional 
constant of the motion. Indeed, from the form of the Euler vector field (j4.3p . we easily ascertain 
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that {TCi, Ci}Pi =0 and a direct computation (or a careful examination of the generahsed Lenard 
relations associated with P, Q) shows that (^i commutes with 7i2 as well. 

So, there must be a functional relation between TCq, C2, Tli, 'H2 and Ci = ~ ^2 alone, that 
is, we expect $1 to be independent of ^1. Taking into account that the elements H are quadratic 
functions of zi , Z2, we look for a relation of the form 

= aCl + Hi + I5n2 + 71^0 + 72C2, 

for some unknown constants a, 7^, i = 1,2. Indeed such a relation can be found with, 
respectively, 

a = 2 ■ , /? = ■ , 7i=/ii+//2, 72 = 0. 

/^i + /^2 /^i + ^2 

To find the second separation relation is slightly more involved; still the idea is to look for a 
relation quadratic in and affine in TIq, C2, Ti-i, 'H2, with coefficients that may depend non 
trivially on the co-ordinate A2, i.e. a relation of the form 

^2 = P{X2)CI + gi(A2)Wl + q2{X2)n2 - ^(A2, Wo, C2), (4.9) 

for some functions p, gi, q2 that depend only on A2, and for an unknown function ^{\2;Ti.o, C2), 
affine in the Cj's. After a direct computation one sees that the problem can be solved, and that 
the second separation relation has the form ()4.9p . with 

gi = A2, (?2 = 1, p = -2/ii(F(u)2G(u)2), ^ = Xlno-fi3F{u)G{u)C2, (4.10) 

where the functions G(u) and -F(u) are given respectively by (j4.6p and ()4.7p . 

It can be noticed that the relation (I4.10p is quadratic in the momentum ^2, and is algebraic - 
rather that polynomial - in the co-ordinate A2, owing to the relation (j4.7p . 

The question whether these techniques might be useful in the study of the general so (4) Euler- 
Manakov top, that is, the case /U3 7^ /i4, is still under investigation. What is still lacking, in this 
general case, is the determination of the suitable deformation vector field Z; it is conceivable 
that a careful use of the theory of elliptic function may provide the answer. 

Also, the connection of the results herewith presented with the Clebsh model, according to 
an isomorphism described in [5] might be worth of further investigation^, as well as the link 
with the setting, within the bihamiltonian theory, of the Separation of Variables problem for 
Lagrange tops presented in |24j . 
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